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Abstract 

A review of the random magnetic impurity model, introduced in 
the context of the Quantum Hall effect, is presented. It models an 
electron moving in a plane and coupled to random Aharonov-Bohm 
vortices carrying a fraction of the flux quantum. Recent results on its 
perturbative expansion are given. In particular, some funny families 
of integrals show up to be related to the Riemann C(3) and C(2). 

1 Introduction : 

The single Aharonov Bohm (A-B) flux line (infinitely thin and im- 
penetrable vortex) piercing the plane at a given point was considered 
in a seminal paper in 1959 Q]. A more complex system consists of 
several A-B lines piercing the plane at different points [2]. A sim- 
plification arises when the locations of the punctures are random : 
consider the Random Magnetic Impurity model [Hj introduced in the 
context of the Integer Quantum Hall Effect. It consists of an electron 
coupled to a Poissonian distribution of A-B lines (the so-called mag- 
netic impurities) having a mean density p and carrying a fraction a 
of the electron flux quantum <E> . Periodicity (a <E [0, 1]) and sym- 
metry considerations (no priviledged orientation to the plane) allow 
to take a e [0,1/2]. The Poissonian disorder defines the infinitesimal 
probability dP(N, f [, f 2 , f ' N ) of finding N impurities at position 
f 2) •••! r n 

i—l 

where p is the mean impurity density : in the thermodynamic limit 
N, V — > oo (V the area of the plane), p =< N > /V. 
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The model has two parameters, a and p. Periodicity a — > a + 1 
and symmetry around a = 1/2 in the spectrum imply that the density 
of states and the partition function are invariant under a — > 1 — a, 
thus depend in fact on p and a(l — a). 

It has been shown via path integral random walk simulations that, 
when a — > 0, the average density of states of the electron narrows 
down to the Landau density of states for the average magnetic field 
< B >= pa^ with broadened Landau levels (weak disorder). On the 
contrary, when a — > 1/2, the density of states has no Landau level 
oscillations and rather exhibits a Lifschitz tail at the bottom of the 
spectrum (strong disorder). 

In the path integral formulation, the average partition function 
rewrites as an average over C, the set of closed Brownian curves of a 
given length t (the inverse temperature) 

<Z>=Z <e"^ s ^ an - 1 ) > {c} 

where S n is the arithmetic area of the n-winding sector (i.e. the sector 
wounded around n times by the path) of a given path in {C} (Z a is 
the free partition function). It amounts to say that the Poissonian 
A-B lines couple to the S n 's, a different (intermediate) situation from 
the single A-B line, which couples to the angle spanned by the path 
around it, and from the homogeneous magnetic field, which couples to 
the algebraic area enclosed by the path. The average partition function 
rewrites as 

<Z>=Z J e- pt< - s+lA 1p(S,A)dSdA 

where S = | ^2 n S n sin 2 (7ran) and A = i S n sin(27ran) are ran- 
dom Brownian loop variables and P(S, A) is their joint probability 
distribution. Since 0] S n scales like t -in fact < S n >= t/(2irn 2 ), even 
more, for n sufficiently large, n 2 S n n 2 S n >= t/(2n)~, the variables 
S and A are indeed t independent. 

In fact, very little is known on the joint distribution function P(S, A): 

i) when a — > 0, < Z >^ Z Q < e i<B> ^^ nSn >{c}, which is,as 
expected, the partition function for the homogeneous mean magnetic 
field < B >, since J2nS n is indeed the algebraic area enclosed by the 
path. 

ii) when a — > 1/2, < Z >^ Z a < e^^odd Sn >{c}, implying 
that, for the average density of states (in terms of p {E) the free density 
of states) 

<p(E) >=p (E) [ P P(S')dS' 
Jo 

Here P(S') is the probability distribution for the random variable S' = 

lJ2n S n, n odd- 
Turning now to a quantum mechanical formulation, the average 
partition function is 
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<Z>=e-^^^f<Z N > 

N 

where < Zjq > is the average N impurity partition function 

r N a~i 
<Z N >= jH^ATre-^ 



for the TV-impurity quantum Hamiltonian H^. 

The partition function of the mean magnetic field (B) = pa<p is 
reproduced as a power series in (pa) n (mean field approximation where 
the local magnetic field B(r) = 4>~Y^ =1 5{r — r' i ) is replaced by its mean 
value). Perturbative corrections p n a m ,m > n to the mean magnetic 
field expansion originate from disorder: in the 1-impurity case the 
correction is known < Z\ >= Z a + a(a — l)/2; in the 2-impurity 
case one encounters non trivial Feynman diagrams at order p 2 a A , i.e. 
an electron interacting with 2 impurities 4 times, at order p 2 a 6 , i.e. an 
electron interacting with 2 impurities 6 times, etc... 

In the symmetric gauge and in configuration space 

\ i=i v l > / i=i 
The nonunitary wavefunction redefinition 

JV 



4>{?) = \{\r f -?rm 

i=l 

results in the Hamiltonian Hjq acting on ip (z complex coordinates) 

N 1 

H N = -2d- z d z - 2a y ^^jd z =H + V N (z) 



1 Z - Z l 
2 — 1 ' 

where quadratic interactions in the vector potential have disappeared. 

The perturbative expansion for the thermal propagator Gp = er^ Htq 
at order n is 

Gp = (-1)" df3 1 d/3 2 --- d(3 n 
Jo Jo Jo 

G p _ Pl V N U 0i _p 2 ■ ■ ■ V N Lrp n 

where G^ = e~ l3H ° is the free propagator. The partition function at 
order n follows as 
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Z$> = TrGf = Tr&f = (-l)"Tr / dpi / dp 2 ■ ■ ■ / df3 n 

Jo Jo Jo 



t> r'(°^ if r<( ) 



In the configuration space one has 



Z ( n ] = (-1)" f dPi f 1 dfc ■ ■ ■ 1 dp n [ d 2 Zl ■■■ [ d 2 z n 
Jo Jo Jo J J 

Vn(zi)G%_p 2 (zi, z 2 ) ■ ■ ■ V N (z n )G { plp n ^ 0i (z n , zx) 



G^(z 1 ,z 2 ) = (z 1 \e-^\z 2 ) = ^e- i ^ 



N 1 

V N ( Zl ) = ( Zl \V N \ Zl ) = -2aJ2^^jd Zl 

Z\ z- 



i=l 



with the average over disorder 



d 2 z'j- — = TTZl (1) 

Zx - z[ 



,2 / 1 / z l , z 2 



d 2 z\ _ _,_ = 7r (^^ + ^±^) (2) 

Zi - z z 2 - z' z 2 - Zl Zl - Z 2 



Note that if at first order in a an ambiguity arises since V{zi) is a 
differential operator acting on the constant G^\zi,zi), this ambiguity 
can be lifted by a long distance harmonic well regularization or by other 
means (see for example p])- 

In momentum space one has (p complex coordinate) 

G ( °\p uP2 ) = (pi\e-P H °\p 2 )=e-^6 2 (pi-p 2 ) 



N 

Pi -P2 



[pi\V N \p 2 ) =-2aV^^e^i-^+' l - c - 

* * Tli Tin 



; = 1 

with the average over disorder implying momentum conservation. 

With these perturbative Feynman rules at hand, one can proceed 
with the perturbative expansion of the average partition function. The 
first non trivial diagram occurs at order p 2 a 4 . Denoting the tempera- 
ture differences as a = Pi — p 2 , b = f3 2 — /3 3 , c = (3 3 — 04, d = P + P4 — Pi , 
the momentum integration yields an integral over the temperatures 

1 f P JO f 01 JO f 132 J3 f^jof 2 (a + c)(b + d) 

4 



i rp rPi rPi rP3 

diag(p 2 a 4 ) = — / dpi / dp 2 / d/3 3 / dp. 
P Jo Jo Jo Jo 



P abc + bed + cda + dab , 
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Integrating on 04, 3 , etc., leads for the non trivial part to [Hj 

v r « r « r ■» ; _ ^ > + cm 



a&c + ocd + cda + dab 16 

(3) 

where £(3) = 7£(3)/2. One can go a bit further : first rewrite this 
integral as 



JO 



abc + bed + cda + dab 



where d = — (a + b + c) is understood. This is 

dadbdedd , d ( a + c )( b + d ) up -(a + b + c + d)) 
a.b,c,d=o a bc + bed + cda + dab 

or, using obvious symmetries, 

f°° ab 

dadbdedd— — — 6(0 - (a + b + c + d)) 

J a b,c.d=o °" c + oca + cda + dab 

Therefore ©becomes 

dadbdedd— — 6(0-(a+b+c+d)) = 



2 1 + C(3) 



abc + bed + cda + dab 16 

(4) 

Then Laplace transform iQJ (J Q d0e /3s ...) to obtain 

1 -s(a+b+c+d) _ 2 1 + C(3) 



da db dc dd — — i ^ — e 

a,fe,c,(i— 



Fix s = 1 and, in the denominator, exponentiate - + \ + - + ^: 

1 ^-(a+b+c+d)-t(^ + i + i- + i,) _o 1 +C( 3 ) 



dadbdedd / dt-,e-( a+b+c+d ™+i + i+^ =2- 

cd 16 



a^b.c.d— 

Change variable u = 2-s/t and use 

a n - l e- a -i = 2K n (u)(^) n 
where the K v (u)'s are modified Bessel functions [§], to finally obtain 

\du(uK x (u)) 2 K Q (uf = 1+ ^ (3) (5) 

16 

which also been derived in a different context 

It can be generalized jH] : the family of integrals 



p„(0000)= / du u n Kq(u) n>0 
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p„(0011)= / duw n+1 if (V) 2 #i(w) 2 n>2 
Jo 

poo 

p n (llll)= / duu^K^u) 4 n > 4 
Jo 

(n even) 
and 

/■OO 

i„(0001)= / duw n+1 Aro(w) 3 ^i(u) n>l 
Jo 

i n (0111)= / duw' l+1 i^ (w)A'i(u) 3 n>3 
Jo 

(n odd) can be shown to be again related to C(3), i.e. of the form 
©• To show this, integrate by parts (use dK (u)/du = —Ki(u) and 
d(uKi(u))/du = —uKq(u)) to obtain 

p„(0000) = — l_i n+1 (0001) n>0 
n + 2 

p„(0011) = -(i„+i(0001) + i„ +1 (0111)) n > 2 
n 

p„(llll) = i„+i(0111) n>4 

n — 2 

which implies 

2np„(0011) = (n + 2)p„(0000) + (n - 2)p n (llll) n > 4 

Therefore z„+i(0001), i„+i(0111), andp„(0011) are obtained fromp„(0000) 
and p n (llll) which remain to be determined: further integration by 
parts gives 

i n (0111) = — !—((p n+1 (llll) + 3p„+i(0011)) n > 3 
n — 1 

i n (0001) = — |— ((p n +i(0000) + 3p„+i(0011)) n > 1 
71 + 1 

thus the recurrence relation acting on a 2 dimensionnal vector space 

p n+2 (0000) \ _ 1 / (n + 2) 2 (5n + 4) -3n 2 (n - 2) \ / p„(0000) 



p„ +2 (llll) y 2 5 (n + 2) V -3(n + 2) 2 (n + 4) n(n - 2)(5n + 16) J \ p„(llll) 

(6) 

for n > 4. For < n < 4 direct computations (one has to rewrite these 
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simple integrals in terms of multiple integrals) give 

C(3) 



Po(0000) 
ii(0001) 
p 2 (0000) 



2 2 
2 3 

-3 + C(3) 
2 4 



p 2 (0011) = 1 + known =© 



* 3 (0001) 



-3 + C(3) 



2 4 

13(0111) = —~ obvious 
I 2 - 

and the n = 4 initial conditions for the recurrence are 

-3 3 + 7C(3) 



p 4 (0000) 



2 6 



53 - 3 2 C(3) 

P4(11H) = (7) 

For example, in terms of the corresponding multiple integral, po(0000) = 
C(3)/2 2 is nothing but 

poo 1 

/ dadbdcdd- — e -(a+b+c+d) = 7 ^gx 

Ja,b,c,d=0 a " C + " C " + C " a + 

If one sets n — 2k and defines (j'fc(O) = p2fc(0000)/(2fc)! and gfe(l) = 
P2/c(lHl)/(2fc)! the recurrence @ becomes 

9fc+ i(0) \ = 1 / (5fc + 2) -3fc 2 (fc - l)/(fc + l) 2 \( q k (0) 
9fc+i(l) J 2 4 (2fc + l) v, -3(fc + 2) fc(fc-l)(5fc + 8)/(fc+l) 2 j ^ ?fc (l) 

(8) 

k > 2. In the asymptotic regime fc — + 00 



J_ / 5 -3 
2 5 I -3 5 



diagonalizes as 



V2^-l 1/2^-3 5 J V2 V 1 1 / 2 4 ^ 4 

with eigenvalues 1/4 and 1/16. It is a remarkable fact that the initial 
conditions 0, that is for the q's 

-27 + 7C(3) 53-9C(3) 
92(0) ~ (3)(2«) ' q2{l) - (3)(28) 

are such that the asymptotic behavior of <fc(0) and <ft(l) happens to 
be governed by the smallest eigenvalue, i.e. lim^oo qk(0) — (1 /16) &— 2 
and Um^oog^l) ~ (l/16) fc - 2 . 
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This is not the end of the story: pushing the perturbative expansion 
at order p 2 a 6 [Hj one encounters the integral J °° du I\ (u)K (u) 2 (uKi(u)) , 
where I\{u) is again a modified Bessel function 0. This integral is 
found to be 



duI 1 (u)K (u) 2 (uK 1 (u)) = ^ (9) 
It is not an isolated case: the family of integrals 

/>oo 

p°(000)=/ duu n {I {u)u)K (u) 3 u 2 



p 1 n (000)= duu n h{u)K (u) 3 u 2 
Jo 

/>oo 

p°(001)= / dM"(/„(u)u)K («) 2 («Xi(ii)) 
Jo 

/>oo 

p*(001) = / du u n h{u)Ko{u) 2 {uKx{u)) 
Jo 

/>oo 

p°(011)= / d U u"(/ (u)u)^o(u)(^i(w)) 2 
Jo 

/>oo 

pi(011)= / d«u n /i(u)lfo(u)(utfi(tO) 2 
Jo 

poo 

p°(ih)= / d«« n (j («)«)(«^i(«)) 3 

JO 

/>00 

pi(lll)= / duti n ii(u)(tififi(u)) 3 



(n > 0, n even) can be shown to be also related 1 to £(2) (note that Q 
is pj(001) = C(2)/8). 

Indeed, integration by part (this exercice is left to the interested 
reader, use dI °( u ' = Ii(u) and ^"^f"^ = uIq(u)) gives that 

P°+ 2 (000) \ _ 1 / (n + 4) 2 (57i + 14) 3(n + 2) 2 (n) \ / p°(000) 



ri +2 (Hl) y 25(n + 4) V 3(n + 4) 2 (n + 6) (n + 2)(n)(5n + 26) ; V ^ (111) 

(10) 

and that all other integrals are given in terms of p„(000) and p 1 (111). 
For example 

2 5 (n + 3)(n + 4)p° +2 (011)= + 8(n + 2) 2 (n + 4) 2 p° (Oil) 

- (n + 4) 2 (3n 2 + 15n + 14)p° (000) 
+ ?i(n + 2)(3n 2 + 27n + 58)p, 1 l (lll) 



1 One also has 

/ du(7 (M)M)i«:oM 3 =3<(2)/8 
Jo 

but, surprisingly, J °° dit Ii(u)Ko(u) 3 is not related to C(2). 
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Amusingly enough, the recurrence iflfljl is identical to © up to n — > 
n + 2 and to the off-diagonal global sign, the origin of which being 
easily tractable to the sign difference in the differentiation relations 
between Iq and I\ on the one hand, and Kq and K\ on the other hand. 
Setting as above n — 2k, k > and redifining 



A (000) 



2(fc + l))! 
%v ' (2(k + l))\ 



the recurrence becomes 



9fc+i(°) A 1 f ( 5fc + ? ) 3fc(fc + l) 2 /(fc + 2) 2 \ / g°(0) 

J 2\2k + 3) ^ 3(fc + 3) k(k + l)(5fc + 13)/(* + 2) 2 j ^ 

with again the eigenvalues 1/4 and 1/(16) for the asymptotic matrix. 
The initial conditions are (by direct computation) 

?o» = C(2) 



9o(0) = -g+3C(2) 

where C(2) = 3C(2)/2 6 . 

To conclude a few remarks can be made: 
-in both recurrences, the building blocks, £(2) = 3C(2)/2 6 and £(3) = 
7£(3)/2, are not pure £(n)'s, but rather (2™ — l)£(n)'s. This indicates 
that one deals with alternated Euler sums, 

p=i ^ 

cm - . ?lzic(«) 

-it has been impossible so far to generalize to £(4), C(5), ■■ in terms of 
simple integrals on the product of modified Bessel functions 2 . May be 
one should look at double nested integrals as indicated j2] by compu- 
tation at order p 2 a 6 . 

2 Note that integrals on a product of four modified Bessel functions seem so far needed: 
for example 

is not simply given in terms of £'s. 



9 



References 



[1] Y. Aharonov and D. Bohm, Phys. Rev. 115 (1959) 485; for an 

earlier description of the same effect see: W. Ehrenberg, R. W. 

Siday, Proc. Phys. Soc. London, sec. B 62 (1949) 8 
[2] P. Stovicek, Phys. Lett. A 142 (1989) 5; S. Mashkevich, J. 

Myrheim and S. Ouvry, Phys. Lett. A 330 (2004) 41 
[3] J. Desbois, C. Furtlehner and S. Ouvry, Nuclear Physics B [FS] 

453 (1995) 759 

[4] Alain Comtet, Jean Desbois and Stephane Ouvry, J. Phys. A 23 
(1990) 3563; Wendelin Werner, These, Universite Paris 7 (1993) 
and Probability Theory and Related Fields (1994) 111 

[5] Alain Comtet, Yvon Georgelin and Stephane Ouvry, J. Phys. A: 
Math. Gen. 22 (1989) 3917 

[6] see Gradshteyn and Ryzhik, Table of Integrals, Series and Prod- 
ucts, Fifth Edition, Academic Press, 8.432 1. and 8.432 6 

[7] L. Samaj and I. Travenec, J. Stat. Phys 101 (2000) 713 

[8] Cyril Furtlehner and Stephane Ouvry, |math-ph/0306004| "Inte- 
grals involving four Macdonald functions and their relation to 
7C(3)/2" 

[9] Cyril Furtlehner and Stephane Ouvry, in preparation 
[10] This was shown in collaboration with Guillaume de Montlaur, 
stage de Maitrise Paris VI (summer 2004) 



10 



